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ABSTRACT 

o _ 

' Context. To determine the local dark matter density around the solar system is a classical problem in astronomy. Recently, Garbari 

et al. have devised a novel method to determine the local dark matter density from stellar distribution and vertical velocity dispersion 
' ^ . profiles perpendicular to the Galactic plane. Their method has the advantages of abolishing conventional approximations and using 

' only a few assumptions. Their determinations, however, preferred relatively high dark matter densities, which are inconsistent with 

other recent works. 

OO Aims. The first aim of this paper is to carefully scrutinize the Garbari et al. method. It should be investigated whether the discrepancy 

(N . between Garbari et al. and the others can be attributed to systematic biases inherent in their method. The second aim is to examine 

influences by observational uncertainties. We discuss how the determinations of the local dark matter density vary with observational 
precisions on parallax, proper motion and line-of-sight velocity measurements. 

Methods. To examine the influences by the observational imprecision, we create mock observation data for stars being dynamical 
tracers based on an analytical galaxy model and apply parametrized observational errors to the mock data. We evaluate the accuracy 
of determining the dark matter density by applying the Garbari et al. method to the mock data. In addition, we estimate a sample size 
,.1^ ' and observational precision to determine the dark matter density with accuracy. 

^L{, Results. We find that the Garbari et al. method is capable of determining the local dark matter density with accuracy if the sample size 

and observational precisions are satisfactory. The required sample size is approximately 6,000 stars. However, the parallax and proper 
motion errors can cause overestimation of the dark matter density. We estimate the required precisions of the parallax measurements 
to be approximately 0. 1-0.3 milliarcseconds; the proper motion precisions do not seem to be as important as the parallax eiTors. Also, 
^ , we find that the line-of-sight velocity errors can cause either underestimation or overestimation of the dark matter density, which is 

contingent on distance-dependence of the errors. 

Conclusions. From these results, we infer that Garbari et al. would be overestimating the local dark matter density due to use of 
imprecise data from the Hipparcos catalog; however, we emphasize the capability of their method. We expect that Gaia will provide 
^ , data precise enough to determine the local dark matter density using the Garbari et al. method. 

\ Key words. Methods: analytical - Astrometry - Galaxy: kinematics and dynamics - solar neighborhood - dark matter 
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. 1. Introduction [T99TI ISofue eTal]l2009t Ifagang et al.ll20Tl : ISofuel (l20T2lP l has 

' , ,. , , . ,. , , recently estimated that pSL = 0.00612 ±0.00080 Mopc-^; 
In recent years, dedicated searches tor candidates or dark mat- , , , , , m _^r7i 
m ter (DM) particles have intensified. Experiments aiming at di- *e other recent works also proposed p'^ < 0.01 M© PC-'B 
^ . rect detection of the DM particles look for signals from recoil of However, the analyses using rotation curves posit g obal niass- 
> ■ DM particles with nuclei inside the detector (e.g., Spooner 2007; I^^^ ^"'^ l^f ^l assumption of the Galactic halo, 
r fBaudis 2012, and references therein). Thus, the eventrates of the Galactic DM halo is far from a round shape, the ro- 
X direct detection are clearly proportional to the local DM density ^^^1°" '^'^''^^ ^fl^'^'^''" '^^'^ !° erroneous determmations 
b (LDMD) around the solar system, p® . However, the event rates ^^"'^^ estimated LDMDs are spherically averaged densities 
- ■ ■ have degeneracy between the LDMD and cross-sections of the (^-S-' l ^empet aL l M Cosmological simulations have m- 
DM particles, therefore we cannot determine the cross-sections deed shown that DM ha los are generally aspherical (e^^ 
without knowledge of the LDMD. This is why the problem of 'Dubinski & CarlberdMBl Allgood et alj2 00i| Bett et al. 2007t 
determining the LDMD at the solar position has recently been ^^^hlen et al. 2007). Moreover, recent studies on dynamical anal- 
atti-acting a great deal of attention. ^ °f *f . Sagittarius stream have preferred a neariy oblate 
Because of the "dark" nature, however, it is difficult to un- ^alo in which the minor axis is approxima tely a igned with 
derstand details of the stmctures and density distribution of the *^ ^^"^ °f ^^g^t to the Galactic center (e.g., | Law & Maiewski| 
Galactic DM halo. Usually, we have to resort to dynamical anal- 
yses of observable phenomena, such as motions of stars, gas and 

satellite galaxies. For example, a rotation velocity curve of the ' ISofuel ( l20l2h used an extrapolation of a Navarro-Frenk- White den- 
Galactic disk is often used to measure the DM halo mass and sity profile, which is determined by fitting to outer Galactic regions, 
density profile (e.g., lAUen & MartosI 119861; lAllen & SantillanI However, he mentioned that his rotation curve inside ~ 20 kpc can- 

not constrain the DM density profile. 

* e-mail: shigeki . inoueOnao .ac.jp ^ I Mq pc^^ ~ 38.0000 GeV cm^-^ 
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20101: llbata etaP 120121: ^eg & Widrowl I2013PI although see 
Oiling & Merrifieldl 120011: iDebattista et al.l l2013h . Thus, more 



accurate determination of the LDMD requires following a more 
direct methodology than modeling the whole halo. 

Some studies have tried to determine the LDMD with anal- 
yses of local kinematics. For example, Salucci et al. (201^ 
used an equation of centrifugal equilibrium and an observed 
local radial gradient of the rotation curve: They determined 



pQ, = 0.01 13 ±0.0055 



Mpj pc . In addition, other methods 



of fitting various dynamical parameters, such as Galactic rota- 
tions, halo masses and thicknesses of gas layers, also derived 
LDMDs of pf^ - 0.01 Mwp c ^ lOlling & Merrifieldj l200d 
1200 It IWeber & de Boe3 l2010t ICatena & UlUo .2010.) . These 
methods relying on the Galactic rotations and its gradients, how- 
ever, still posit axisymmetric assumptions. 

As another approach, vertical motions of disk stars near 
the sun are frequently used to determine the LDMD. The 
Galactic potential of visible and invisible matters (i.e., the to- 
tal density) can be derived from zenithal profiles of density 
and velocity dispersion of tracer s tars via the Boltzmann and 
Poisson (or Jeans) equatio ns (e.g., IOortlll932t lBahcalll[T984bt 
iKuiiken & Gilmorall989cl) . The visible matter density is mea- 
sured by observations and can be subtracted from the total den- 
sity, and then one can determine the invisible matter density. 
This method is advantageous in terms of being independent 
from the global shape of the DM halo. Previous studies have 
determined that < 0.01 Mp pc~^ Bienayme et al. 



1987 



1991 



Creze et al. 1989; Kuiiken & Gilmore 1989a b; Kuiiken 

iFlvnn & Fuchs 1994 ;' Pham 1997; Creze et al. 1998; 

Holmberg & Flvnn 2004; Bow & Tre maine 2012; Zha ng et al.l 
2OL20 although some earl ier works indicated a prefer ence for 
much higher LDMDs (e.g.. lOort|[T932. 1960; Hill 1960; Bahcalll 

[T984aiicl) . 

Recently, iGarbari et all (l201l[ |2012[ hereafter Gil and 
G12, respectively) have devised a new method to determine 
the LDMD from stellar motions perpendicular to the Galactic 
disk. Their method was named the minimal assumption (MA) 
method and uses, as the name suggests, only a few assump- 
tions and a flexible 15-component galaxy model (see Sect. |2]). 
With the MA method, they determined the LDMDs of p'^^ — 



hO.008 
-0.009 



M0 pc-^^ in Gil and p^ = 0.022;°;JJ|^ Mq pc 



,-3 



0.033 

in G12. Their results are higher than the conventional LDMD 
by a factor of two or more. They suggested that the discrep- 
ancy from the low LDMD derived from Gal actic ro tations can 
be indicative of existence of a "dark disk' ' (iRead et al. 2009). 
However, IZhang et"al] (1201 2|) more recently determined pf^ = 
0.0065 ± 0.0023 M© pc"^ with a similar but different method, 
whose sample stars are much larger in number than those of 
Gil and G12 and carefully selected in an [Fe/H1-fa/Fe1 dis- 
tribution. Since the LDMD determined by IZhang et all (1201 2h 
is comparable to the estimation from rotation curves, they sug- 
gested that the Galactic DM halo would be almost spherical 
and there is no need to consid e r the existence of a dark disk. 
In addition, iBovv & Tremaind (1201 2h also determined p^ = 
0.008 ± 0.003 M0 pc"^ from kinematics of thick disk stars. 



^ Some studies are, however, inconsistent with the r esult of the oblate 
halo, preferring an approximately spherical halo (e.g.. Ibat a et alj200ll : 
[Pellhauer et al. 2006) and a prolate halo (e.g., Helmi 2004j). 

^ [Moni Bidi n et alj 120121) proposed absence of DM around the sun 
from their analysis of thick disk stars. However, iBovv & Tremaind 
( 120 12f) claimed that this result significantly underestimates the LDMD 
because of a false assumption. 



The MA method of iGarbari et al.l relying on only a few 
assumptions seems to be a promising method, although their 
results are discrepant from th e conventional LDMD. T hey 
seemed to show tha t metho ds of iHolmberg & Flvnnl (l2004l) and 
iKuiiken & GilmorS (Il989alfbll3) have systematic biases to under- 
estimate the LDMD due to unsuitable assumptions of separa- 
ble distribution function and isothermality for tracer stars. The 
MA method, on the other hand, does not need these assump- 
tions. Moreover, they applied the MA method to A'-body sim- 
ulation data, and demonstrated that their method can determine 
the LDMD within the 90 % confidence level. Howev er, as men- 
tioned above, the more recent works bv IZhang et al.l (20 1 2 ) and 
IBovv & Tremaind (12012 ). which also exclude these disputed as- 
sumptions, support the convention al lower LDM D. What is the 
cause of the inconsistency between FGarbari et al .l and the others? 
As our first aim, we examine the capabiUty of the MA method 
by applying it to mock observation data generated with an ana- 
lytical galaxy model. 

Fortunately, launches of nex t-generation astrometr y satel- 
lites are now approaching: Gaia dPerrvman et al.ll2001l) . Nano- 
JASMINE (J apan Astrom etry SateUite Mission for INfrared 
Exploration. lGoudall201 ll) and /MAPS' (Joint MiUi- Arcsecond 



Gaume et al ]|2009») . These satellites will pro- 
stars with unprece- 



Pathfinder Survey. 

vide astrometric information of Galactic 



dented precisions. Therefore, it is worth performing a feasibil- 
ity study of determining the LDMD with data provided by them. 
We study the dependence of a determined LDMD on a sampling 
region and observational errors such as parallaxes, proper mo- 
tions and line-of-sight velocities (LOSVs). Our second aim is 
to evaluate required sample size and observational precisions to 
determine the LDMD with accuracy. 

The basic concept of the MA method is described in Sect. |2] 
Our galaxy model for the mock observations and how to create 
the mock data are explained in Sect. [3] and |4] respectively. A 
Markov chain Monte Carlo (MCMC) technique combined with 
the MA method is explained in Sect. |5] Our results are shown in 
Sect. |6] Discussion and our conclusions are presented in Sect. [T] 
and Sect.|8] respectively. Since Sect.|2]and|5]basically follow the 
same course of action as G12, readers who are already familiar 
with the MA method or only intended in our results may prefer 
to skip these sections. 

2. The MA method 

Here we explain procedures to derive the LDMD with the MA 
method; however, we recommend that readers refer to Gl 1 and 
G12 for the original idea. 

2.1. Basic equations and deduced gaiaxy modeis for fhe MA 
metliod 

To begin with, we assume that the Galaxy can be modeled as 
a superposition of multiple visible components, such as stel- 
lar and gaseous disks, stellar halo, and a dark component, i.e., 
Ptot = L Pi + Pdm, where p,- represents mass density of the i-th 
visible component, and ptot and p^^ are the total and DM mass 
density, respectively. Next, the MA method introduces only three 
assumptions: 

(1) the system is in equilibrium, 

(2) a "tilt" term in the Jeans equation is negligible, 

(3) the DM density is constant in a region we consider 

Gil have confirmed that all of these assumptions hold in an A^- 
body simulation of a barred spiral galaxy. 
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From the assumption (1), the i-th component satisfies the 
Jeans equation. Considering the direction perpendicular to the 
Galactic disk (z-axis), the Jeans equation is described as 



0, 



(1) 



where V,- is density of the ;-th component under the total gravita- 
tional potential, <I>, and the over-line means the averaged value, 
and (7- is vertical velocity dispersion. The first term is often 
called the tilt term. From the assumption (2), Eq. ^ ignoring 
the tilt term can be reduced to 



dVi 




(2) 



Integrating this equation, one can obtain a density profile of the 
i-th population as 



v,-(0) aliz) 



exp 



1 d4> 

al-{z') d? 



dz' 



(3) 



(Jz.,(0) are given as parameters, whereas the profiles of C7V,,(z) 
must be assumed to solve this equation^ Gil introduced 
parametrized runs of C5;,, (z), and G12 assumed isothermality for 
all components of the galaxy model. Following G12, we also 
assume that the model components are isothermal; however we 
discuss the influence of different <y^j{z) profiles on the calcu- 
lation in Appendix lAl With the isothermal assumption, Eq. (|3]l 
becomes 



(7.2.(0) 

^P;(0)-3-7Texp 



(4) 



Here, p,- is mass density: p/ = V,m,-, where m,- is the mass-to- 
light ratio of the i-th population. p,(0) are given as parame- 
ters. Finally, one can describe the profile of the total observ- 
able density as Ps(z) — LiPi(z)- Then, the total matter density 
is Ptot(z) = ps{z) + Pdm from the assumption (3). In addition, the 
total surface density of the observable matters is derived as 



Es(z)=2 / 'ps{z')dz'. 
Jo 



(5) 



Next, ptot is connected to the total gravitational potential, 4>, 
via the Poisson equation: 



d^-^ 1 dV;r{R) 



~ dz^^R dR 
This is transformed to 



dz^ 



47tG[ Ps(z)+Pdn 



1 dv^jR) 

A%GR dR 



(6) 



(7) 



This equation implies that the radial gradient of circular veloc- 
ity, causes a pseudo-mass which can be either positive or 
negative. However, if the z-dependence of this term is ignored, 
it can be simply subtracted from the "effective" LDMD, Peff = 
1 dvl{R)_^ fm?i\\y. Although iBovv & Treniairig (l20T2 h 



Pdi 



AkGR dR 



pointed out that neglecting the z-dependence of can lead to 



^ It is impossible to observe OV,, (z) of all components up to high z, 
therefore we have to assume their runs. We consider the situation where 
only a velocity dispersion profile of a tracer population is available. 



underestimation of the LDMD, they also evaluated the underesti- 
mation as being only a tiny fraction of the matter density at low 
z where the thin disk components are dominant. Additionally, 
this term is observationally accessible via the Oort constants^ 
Accordingly, we exclude the term of circular velocity gradient: 
Eq. O simply becomes 



:4;rG(ps(z) + Pdm)- 



(8) 



One can now compute the potential of the model galaxy, <I>, 
by numerical integration of Eq. (|4]i and (|8) with a given param- 
eter set of p, (0), C7V,,(0) and pdm and boundary conditions of 
<I>=^=Oatz = 0. The total number of the parameters is, 
therefore, 2n + 1, where n is the number of visible components 
in the galaxy model. In this paper, we refer to the computed po- 
tential, «I>, as the "deduced model". 



2.2. A tracer population 

The reliability of a trial parameter set used in the above compu- 
tation of <I> must be evaluated by comparing with observations. 
To this end, the MA method prepares an observed sample of 
"tracer" population stars. The tracer stars must share the same 
kinematic state. The prefeiTed tracer stars are considered to be 
old (well-mixed), bright (observationally precise and accurate) 
and large in number (statistically reliable). 

Let us now assume that a sample of tracer stars is available, 
and we obtain density and velocity dispersion profiles, V(°2^(z) 
and (J°|iJ.|(, (z), in a range from Zmin to Zmax- Substituting (7°]?rlc (z) 
for a, J in Eq. Q, one can deduce a fall-off of the tracer density 
profile under the modeled potential, <I>, as 



wded 



iz) 



fJl obs 
*-'z,trac 



(^min ) 



^trac (^min ) 



^2 obs 



exp - 



1 



dO 



Jl obs 



(zO dz' 



dz' . (9) 



This deduced density fall-off can be compared with the ob- 
served fall-off, Vt°ac(z)/Vtracfe-nin). ^nd goodness-of-fit between 
them is calculated to estimate the values of the model param- 
eters. Although our study uses a single tracer population, the 
MA method does not necessarily require only a single tracer, but 
multiple tracers can also be used. 

Since the MA method usually introduces a large number of 
parameters, an MCMC technique should be applied to explore 
the vast parameter space (see Sect. |5]l. Then, as a result, proba- 
bility distribution functions (PDFs) and their degrees of uncer- 
tainties are obtained for all of the model parameters. 



3. An assumed Galaxy model for mock 
observations 

In this paper, we do not make use of observational data or intend 
to determine the LDMD in the real Galaxy. Instead, our study 
performs mock observations of a tracer sample using an analyt- 
ical model, and aims to assess the capability of the MA method 
by deriving an LDMD assumed arbitrarily. To this end, we pre- 
pare an assumed galaxy model for the purpose of creating the 
mock tracer data. 



* iBovv et all l l2012ah estimated that the contribution of this term to 
the local density is only 0.0002+y Qyjs pc"^, and Gil estimated 
-0.0033 ±0.00 50 Mq pc^^ from observations of iMignardI fcOOOl) and 
lBranhanJ(l2010h . 
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First, we assume that density profiles of all visible compo- 
nents are represented as pf^^ °^ sech^ (z/(2/!,)) although this may 
not be suitable for gas or stellar halo components. Then, the total 
density structure perpendicular to the disk plane is 



(10) 



P.otW=Lpr(0)sech2(^)+Pd„ 



where pf'*'*(0), /i, and are set arbitrarily. Pd^ is a constant ac- 
cording to the assumption (2). By integrating the Poisson equa- 



tion with the boundary condition of ^*fl = at z = 0, the dif- 
ferential of this galactic potential is calculated analytically: 



d*^ 



dz 



-A%G 



^(O)tanh(- 



- PdmZ 



This is related to the total surface density as — A-nG'L^^{z), 
and surface density of all observable matters is Ef'*(z) = 
2^,/!,pf^^(0)tanh (z/(2/i,)). Next, we may integrate Eq. ©from 
arbitrary z to z = +oo. Since pf^^ (z) as z ^ +°°, 



(11) 



pr 



-d4>^ 



dz' 



-dz'. 



(12) 



Substituting Eq. ( fTTT ) and z = in this equation, we can compute 
Oi,r(0) when a set of pf'*^(0), /i; and pdm are given. We refer 
to this galaxy model as an "assumed model". In this study, we 
presume that this assumed model represents the real Galaxy. 

Following Garbaii et al., we assume that the model consists 
of fifteen visible components wi th pr'^(O) hsted i n Table [1] that 
are taken from observations of Flvnn et al.l (l2006l) . The LDMD 
is set to Pdm = 0.01 Mq pc^^. In addition, we determine hi as 
follows: In the calculation above, OV,,(0) are the output, and hi 
are the input parameters; however, cr;,,(0) can be determined 
by observations of nearby stars and gas, whereas /i, are poorly 
known. Therefore, it is preferable to search for a set of /i, which 
matches C7,.,(0) computed by Eq. (fTST i to the observed values. 
Accordingly, we employ an MCMC method to find such a set of 
hi. The result of the set of /i, is shown in the last column of Table 
[TQ After hi are obtained, we recompute o;,, (0) by Eq. (fT2l i with 
the set of hi. These are shown in the fifth column, indicating ex- 
cellent agreement with the observations (the fourth column). We 
adopt pf^^(O) in the third column and CT;., (0) in the fifth column 
for the assumed model described by Eq. (fTOt . 

The total observable surface density of the assumed model 
is Ef'*|,=oo — 49.6 M0 pc^^, this is in excellent agreement 
with the ob served value of E"''^ = 49.4 ± 4.6 Mq pc"^ by 
iFlvnn et al] (^006); however, it should be noted that hi seem to 
be somewhat shorter than observed sca le heights of :^ 300 pc 
for Galactic thin disk components (e.g., Gilmore & Reidlll983t 



iBinnev & Merrifieldlll998t[Bovv et al.ll2012blldl) . 



4. The mock data 

As mentioned in Sect. 12.21 the MA method needs density fall- 
off and velocity dispersion profiles of an observed tracer popu- 
lation. Accordingly, we create tracer stars that are dynamically 
consistent with the assumed galaxy model in Sect. [3] We assign 



^ This MCMC calculation can also be applied to the case in which 
Pdm is taken as another parameter. In this case, the MCMC can esti- 
mate Pdm directly. However, the result prefers an extremely low LDMD, 
Pdm — 0, in such a case. This may imply that the density profiles of 
sech^ (z/(2/i,)) cannot be applied to all components in the Galaxy. 



three-dimensional positions and velocities in the Caitesian co- 
ordinates, {x,y,z,Vx,Vy,u^), to the tracer stars. In this model, the 
sun is presumed to be located at {x,y,z) — (0,0,0) and at rest. 

As well as the model components, we assume the tracer den- 
sity profile to be 



Vti-ac(z) = Vtrac(0)sech2 



2/!„ 



(13) 



where Vtrac(O) is a constant, and we set /itrac — 200 pc. We as- 
sume that the tracer stars are uniformly distributed on xy-plane 
at a height above the plane: Vnac depends only on z. Thus, we 
can place tracer stars at random in spatial coordinates according 
to Eq. ( fT3] ). 

The vertical velocity dispersion profile can be calculated by 
Eq. (fT2l) using the assumed model: 



1 r d<P''\ , 

^{z) = / Vtrac— rrdz • 

Vtraciz dz' 



(14) 



Motions of the tracers in the Galactic radial (x-) and azimuthal 
(y-) dir ections can be take n from the observational fitting func- 
tions of iBond et all (l2010h . We assume that there are no merid- 
ional motions: y^.trac = 0, and velocity dispersion in the radial 
direction is 



1.5 



km s 



(15) 



Mean azimuthal velocity with respect to the sun decreases with 
increasing z. 



V. 



7,trac 



■(z) = -19.2(^)"kms 



-1 



(16) 



This is often r eferred to as asymmetric drift (e.g., 

iBinnev & Merrifield 1998 : Binnev & Tr emain aifOOSl) . Also, the 
azimuthal velocity dispersion profile is given as 



CTy,trac(z) -30 + 3.0 



\ 2.0 



(17) 



Assuming Gaussian velocity distributions with the dispersions 
above, we can assign a velocity vector to each tracer star. 
However, it should be noted that velocity distributions may not 
be Gaussian in the real Galaxy. In addition, the observations of 
Bond et al. (2010,,) are for blue stars, whereas the tracer popula- 
tion should be old (red) stars. 

Now that we can assign {x,y,z,Vx,Vy,v^) to the tracer stars 
according to the equations and assumptions above, the positions 
and the velocities are converted to the spherical coordinates of 
{d,9,(l),vios,fg,u^) which are centered at the sun, where d is 
distance from the sun, and i>ios is LOSV. We apply mock ob- 
servational errors to the stellar positions and velocities (see Sect. 
|6]l, then we turn the errored positions and velocities back to the 
Cartesian coordinates (x',y,z',y'^, f^,, u'). Next, we pick out stars 
contained in a certain sampling region, then we divide the sam- 
ple stars into z-bins; the number of the bins is set to Wbin = 10 in 
this study. The z-bins have the equal widths of (zmax — -Zmin ) / "bin, 
which are centered at Zj. Finally, we obtain profiles of V°^l{zj) 
and C7°]iJ!i(, (zy ) affected by the mock errors. We apply the MA 
method to these errored profiles and examine the influence of 
the observational errors, sample sizes and sampling regions on 
determination of the LDMD. We hereafter set Zmin — 2 kpc and 
Zmax = 1-2 kpc unless otherwise stated. 
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Table 1. Parameters used in the assumed gal axy model. Th is model assumes the fifteen visible compo- 
nents shown in the second column, followin dGarbari et al . in which their disk mass model is taken from 
[Flvnn et al. (2006, the third and fourth columns). In this model, we adopt p^m = 0.01 Mq pc^^ and de- 
rive hi (the last column) with MCMC calculation. Using the set of /?,, the fifth column is recalculated by 
Eq. (fT2l l. Our assumed model of Eg. (T TOll uses the third and sixth columns. The uncertainties on OV , (0) 
are set to be the same as iGarbari et al.l 



# 


Component 


Pi (0) IM0 pc ] 
(observation) 


0^,i(u) [kms I 
(observation) 


OV,,(Uj [Kms J 
(MLML) 


hi [pcj 
(MCMC) 


1 


H2 


0.021 


4.0 


3.80±1.0 


32.4 


2 


HI(1) 


0.016 


7.0 


7.00±1.0 


63.4 


3 


HI(2) 


0.012 


9.0 


9.02±1.0 


84.7 


4 


Warm gas 


0.0009 


40.0 


40.1 ±2.0 


564 


5 


Giants 


0.0006 


20.0 


20.0±2.0 


224 


6 


M„ < 2.5 


0.0031 


7.5 


7.56±2.0 


69.2 


7 


2.5 <M„ < 3.0 


0.0015 


10.5 


10.7±2.0 


103 


8 


3.0 < M„ < 4.0 


0.0020 


14.0 


14.1±2.0 


144 


9 


4.0 < M„ < 5.0 


0.0022 


18.0 


18.1±2.0 


197 


10 


5.0 <M„ < 8.0 


0.007 


18.5 


18.5±2.0 


203 


11 


M„ > 8.0 


0.0135 


18.5 


18.4±2.0 


201 


12 


White dwarfs 


0.006 


20.0 


20.1 ±5.0 


226 


13 


Brown dwarfs 


0.002 


20.0 


20.0±5.0 


224 


14 


Thick disk 


0.0035 


37.0 


37.7±5.0 


519 


15 


Stellar halo 


0.0001 


100.0 


99.8±10.0 


1891 



5. The MCMC method 

The MA method is combined with the MCMC technique to ex- 
plore the vast parameter space. We introduce fitting errors of the 
tracer density bins as 



^trac i^j)/ ^trac (^min ) Vtac i^j) / ^ttac (Zmin ) 



ev, 



obs 



fe)/^trac(^mm) 



,(18) 



and an eiTor of the suiface density of all observable matters as 

2 



^.v (-^max ) ^.v (-^max ) 



£^.v (-^max) 



(19) 



where e is a constant that mimics uncertainty in actual observa- 
tions (e.g., stellar mass estimation eiTor). Throughout this study, 
we set e = 0.1. Although the setting of e is arbitrary, gener- 
ally the result does not depend too much on e. Ef^(zniax) is 
47.1 M0 pc^^ in the assumed model. We define the total error 
of the fittings to be 



Xtot ' 



(20) 



where w is a weight of the suiface density eiTor; Q < w < 1. 
A larger w means that the surface density fitting has stronger 
constraint in determining the best-fit parameters. 

We use a Gaussian function as our likelihood function: P = 
sxp(— ziot)- ^ MCMC chain computes the likelihood at each 
step and decides whether a proposed next step is accepted or 
rejected according to a ratio of the likelihoods (e.g., Gregor^ 
I2OIO ). We use a Metropolis-Hasting algorithm for our sampling 
scheme and run the MCMC chains of 50,000 steps: however, the 
first 5,000 steps are excluded as a bum-in period. 

Parameter ranges surveyed in our study are from to 
0.2 M0 pc-^ for pdm, ±10 % and ±25 % for pf '(0) of stel- 
lar and gas components, respectively. The ranges for C7,,,(0) are 
listed in Table [T] and are much wider than actual observational 



uncertainties (iHolmberg & Flvnnll2000l) . The medians of resul- 
tant PDFs are used as the best fit values. The 90 % confidence 
levels (intervals between the 5th and 95th percentiles) are used 
as ranges of uncertainties of the calculation. 

6. Results 

6.1. Application to tiie analytical solutions 

Before adapting the MA method to the mock observation data, 
we directly substitute the analytical velocity dispersion profile 
of Eq. ([Ill in Eq. ^ and run the MCMC by fitting with the 
analytical density of Eq. (T3[ . In such a case, we can investigate 
intrinsic systematic biases of the MA method. 

To begin with, we study the case of w = 0: no constraint by 
the surface density. The result is shown in Figure [T] We can see 
that the LDMD is accurately determined; besides, the local and 
suiface observable densities obtained are accurate enough to be 
within the 90 % confidence levels. From this result, we expect 
that the isothermality assumed in Eq. (|3]l for the deduced model 
does not lead to significant systematic biases (see also Appendix 
lAb and that the MA method would be capable of determining the 
LDMD accurately if there are no observational errors and a sam- 
ple size is large enough. We find, however, that p, (0) and av,, (0) 
of each component are largely degenerate with one another and 
not particularly well constrained. 

6.1 .1 . The surface density constraint in tlie MCIVIC 

In the right panel of the Figure[T] the shape of the PDF is inclined 
from the upper left to the lower right. This implies that pdm and 
Es have a correlation. Therefore, it may be expected that if Eg can 
be determined more precisely, the error range of Pd^ can become 
narrower: the constraint on Es in the MCMC computation may 
help to suppress the uncertainty of p^m. 

Figure |2] shows our results in the cases of w = 0.1, 0.5 and 
0.7. Contrary to the expectation above, the error ranges of pdm 
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Fig. 1. PDFs of LDMD, local and surface densities of the all vis- 
ible matters. Here, we adopt the analytical solutions of the tracer 
density and velocity dispersion profiles and w = in the calcula- 
tion. The left and right panels illustrate the PDFs of LDMD v.s. 
local and surface densities of the visible matters, respectively. 
The brighter regions mean higher probabilities. The green solid 
lines indicate the ranges of the 90 % confidence levels, and their 
intersection points are the medians. The red dotted lines indicate 
the true values in the assumed model. 
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Fig. 2. Same as in Figure [T] but adopting w 
from top to bottom. 
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Fig. 3. Sampling regions for the tracer stars. These are axisym- 
metric. 



widen as w increases; in addition, the inclinations of the PDFs 
in the right panels become steep. This is due to the constraint 
on the tracer density fall-off becoming weaker as w becomes 
larger; therefore, too large w should not be used. The best-fit 
values of Pdm, and Ps (0) are little affected by w. The PDFs 
in the case of w = 0.1 differ very little from those in the case 
of vv = (Figure [U although the error range of Ps(0) is some- 
what narrower Hence, we consider that the constraint on Eg in 
the MCMC may not necessarily make the determination better 
Following G12, we hereafter set w = 0.10 

6.2. Application to tlie mock observation data 

We generate the mock tracer population following the proce- 
dure described in Sect. |4] and observe the sample stars con- 
tained in various regions. We test five sampling regions: "nar- 
row cone", "very narrow cone", "wide cone", "narrow cylinder" 
and "wide cylinder". These cone regions are directed to the z- 
axis (0=0) and contain stars within 9 < 22° (narrow), 2.4° 
(very narrow) and 40° (wide). The cylinders cover the regions of 
R = vx2+y2 Q 5 ]jp(; (narrow) and < 1 .0 kpc (wide). Figure 
|3]dehneates these sampling regions. 

6.2.1. Sample sizes 

If tracer stars are insufficient in their sample size, statistical un- 
certainty may lead to an ill-determined LDMD. Here we investi- 
gate the sample size required to determine the LDMD precisely. 
For this reason, we assume no observational errors here. In this 
case, since our tracer model assumes the horizontally uniform 
distribution and kinematics, the widths of the sampling regions 
do not change the results. The cylinder sampling region con- 
tains more stars at lower z since the stellar density decreases 
with increasing z, whereas the cone region can cover wider areas 
at higher z. Therefore, if the sample size is fixed, samples in the 
cylinder and the cone can be statistically reliable in low and high 
z regions, respectively. 

We perform the calculations using the MA method with 
tracer samples of 1,500-24,000 stars. Our results in the cases of 



^ InG12,w = 0.1 was adopted although they did not introduce the 
parameter explicitly: they used nine bins of a tracer density profile and 
a surface density. 
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0.03 




Fig. 4. Results of derived LDMDs with various sample sizes. For 
each sample size, we performed the same computation six times 
with samples generated by different random seeds. The points 
and the error bars indicate the best-fit values and the 90 % con- 
fidence levels, respectively. The horizontal dotted line indicates 
the true LDMD, Pdm = 0.01 Mq pc"^. 



cone and cylinder regions are shown in Figure For each case, 
we run the same computation six times using different random 
seeds in generating the samples. In both cases of the cone and 
cyUnder sampling regions, the resultant LDMDs with < 3,000 
stars are not precise, and the true LDMDs are often out of the 
error ranges. The samples of > 6,000 stars seem sufficient for 
determining the LDMD within the errors in all runs. Therefore 
we consider that the minimum sample size to precisely deter- 
mine the LDMD would be 6,000 stars. The scatter of the LDMDs 
seems slightly smaller in the case of the cylinder region than in 
the case of the cone region, implying that statistics in lower z re- 
gions are more important. To suppress the statistical uncertainty, 
we set the tracer sample size to be 24,000 stars hereafter, unless 
otherwise stated. 

6.2.2. Distance errors 

Astrometric observations usually measure stellar parallaxes and 
proper motions on the celestial sphere. The parallax measure- 
ments determine trigonometric distances to the stars. Using the 
distances measured, the proper motion measurements can de- 
termine transverse velocities of the stars. One of great advan- 
tages of astrometric obser vations is that systemati c errors are 
generally very small (e .g., Perrvman et al.ll997: de Bruiine et al] 
l2005tlPerrvniaiJ2012l) . For this reason, we examine the influence 
of distance-dependent random errors of parallax measurements. 
We assume no systematic errors — in other words, the mock ob- 
servations are accurate but imprecise — and no errors on stellar 
position measurements (0 and (p); besides, we do not take into 
account dust extinction or binary stars. 

To begin with, we simply analyze how precisions of astro- 
metric measurements varies with distances of stars. Let CJ and 
jj-g be parallax and proper motion in the -direction of a star. 
Apart from bright sources, uncertainty of the parallax is thought 



to be domi nated by photon statisti cs to determine an image cen- 
troid (e.g.. lde Bruiine et al.ll2005l) . In this case, the parallax er- 
ror, £(57, varies with the observed stellar flux, /, as en; 1 / v7- 
Since / °^ d^^, em °^ d. On the other hand, if the source is bright 
enough to obtain a sufficient number of photons and the pho- 
tons do not saturate the image, the parallax error is expected to 
be independent from the distance. Accordingly, we can formu- 
late the parallax measurement error as Ea; °= d", where a = 
or 1. Besides, the parallax itself decreases as distance increases: 
UJ = l/d. This means that the fractional parallax error depends 
on the distance as e^/tiJ (e.g.. Bailer- Jones 2009). Since 

we can equate this fractional parallax error with the fractional 
distance error (FDE), we formalize as 



epDE =A{ - — 
.kpc 



a+l 



(21) 



where A is a parameter which corresponds to the precisions of 
parallax measurements at li = 1 kpc in milliarcsecond (mas) 
units. We apply this FDE with Gaussian probability distribu- 
tions the dispersion of which is Epde to the mock tracer stars. 
In addition, since astrometric measurements determine trans- 
verse velocities as ug — d x jig, transverse velocities are also 
subject to the distance error as v'g = vg^, where v'g and d' are 
errored transverse velocities and distances. However, it should 
be noted that the formulation above may be simplistic, a would 
not be a constant but varies with d in actual observations: when 
a = 1 , the parallax uncertainties cannot be zero atd = Q but are 
limited to calibration errors and/or instrument stability at a cer- 
tain distance. Moreover, in the cases of conical sampling regions 
with fixed angles, observations may suffer from L utz-KeUc er bias 
(Lu tz & Kelkei]ll973HBinnev & Merrifieldll 19981 see also lSmithI 
i2003l) . In this subsection, we examine the influence on determin- 
ing the LDMD by the distance error alone; therefore, we ignore 
the errors on proper motions and LOS Vs here. 

Figure |5] shows our results of LDMDs determined in cases 
of various A, a and the sampling regions. Our results obviously 
demonstrate that the distance errors cause systematic overesti- 
mation of Pdm, regardless of the sampling regions, p^yn is signifi- 
cantly overestimated beyond the error ranges when A > 0.3 in all 
cases except the wide cylinder with a = 1 . Therefore, the paral- 
lax measurement precisions must be A < 0.1 which corresponds 
to < 0.1 mas at c/ = 1 kpc. Only in the case of the wide cylinder 
with a = I, the required accuracy seems to be 0.05 mas. The 
magnitudes of the overestimation seem to depend only weakly 
on a and the shapes of the sampling regions. 

The overestimation shown in Figure |5] is attributed to 
changes in density and velocity dispersion profiles of the tracer 
due to the distance errors. In Figure |6] we show how the dis- 
tance errors deteriorate the reproducibility of these profiles. The 
density fall-offs are similarly affected in both the cone and the 
cylinder sampling regions; however, the velocity dispersion pro- 
files are more significantly changed in the cases of the cylinders. 
In the top panel, the velocity dispersion profile is hardly changed 
in the case of the very narrow cone; however, the LDMD is sig- 
nificantly overestimated even in this case (Figure |5]l. Therefore, 
it can be said that the overestimation can be caused solely by 
the change of density fall-off; however, it can also be caused 
solely by the change of velocity dispersion (see Sect. l6.2.4l i. The 
changes of velocity dispersion profiles are more serious in the 
cases of cylinders than in the cones. This is due to inclusion 
among tracer stars in the cylinder regions stars whose 9 are large 
in low z regions. On such stars, their transverse velocities have 
more information about than their LOSVs. 
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Fig. 6. Changes of density fall-off and vertical velocity disper- 
sion profiles of the tracers. Here we take into account the dis- 
tance errors only. The line colors indicate results of different 
sampling regions and a. The solid and dashed lines indicate re- 
sults of A = 0.3 and 0.5, respectively. Here we used mock data 
of one million tracer stars in all cases. 

6.2.3. Proper motion errors 

Proper motion uncertainties also varies with distance in the same 
way as parallaxes, we therefore formalize the proper motion er- 
ror in the -direction, e^, as 

eM=-Bfj^j masyr-i. (22) 

where B is a parameter that corresponds to the precisions 
of proper motion measurements at li = 1 kpc in the unit of 
masyr^', a is the same as in Eq. (|2TI) . We assume Gaussian 
errors with the dispersion of and add the errors to the true 
proper motions. Errored transverse velocities are v'g = d' x jig, 
where /ig is the errored proper motion. 

In most cases in Figure |5] the LDMD is overestimated by 
the distance errors when A > 0.3; therefore, we here set A = 0. 1 
in all cases. Figure |7] shows our results of adopting various B, a 
and sampling regions. In the figure, the overestimation of LDMD 
is little seen where the cone sampling regions were used, even 
when j3 =3.0. However, slight overestimation can be seen in 
use of the narrow cylindrical region when B > 2.0, and signif- 
icant overestimation occurs in use of the wide cylinder region 
when B > 1.0. Again, this is because the cylindrical regions in- 
clude stars with larger 9 in lower z regions. Therefore, sampling 
in the cylindrical regions seems to have a risk of tracer kinemat- 
ics becoming susceptible to the proper motion errors. Moreover, 
the magnitude of the overestimate seems to strongly depend on 
widths of the cylinders. Accordingly, we suggest that the shapes 
and widths of sampling regions be chosen carefully. The wide 



cylinder would be useful to gather more stars in low z regions; 
however, if precisions of proper motion measurements are in- 
sufficient in actual observations, the tracer sampling should be 
conducted in a conical region. 

In Figure [8] we show comparison of changes of caused 
by the errors on distances and the proper motions. In each three 
panels, the top left and top right panels indicate the influence of 
the errors on distances and proper motions, respectively. If there 
are no errors, the dots align diagonally across the panels. The 
scatter of the dots indicates significance of the error propaga- 
tion to y,. As seen in the figure, the changes of by the proper 
motion errors appear less significant than those by the distance 
errors when A — B. This means that the distance errors are more 
crucial for determining the LDMD than the proper motion er- 
rors; moreover, the distance errors can also change the density 
profiles. The changes of in the very narrow cone are much 
smaller than in the narrow cone, since contribution of the trans- 
verse velocities to u- is ug sin 9. 

From Figure |7]and the discussion above, it may be expected 
that the proper motion measurements are not indispensable to 
derive the LDMD with the MA method if the sampling region 
is a cone. Distances to stars can be measured not only by as- 
trometry but also by other observational means (e.g., photomet- 
ric distances using a color-magnitude diagram). Hence, if the 
LDMD can be determined without proper motion measurements, 
it could be expected that we do not necessarily need astrometry. 
Accordingly, we also examine the cases of ignoring the proper 
motion measurements, where we set /Xe = for all tracer stars. 
Figure |9] shows our results in the cases of the very narrow and 
narrow cone sampling regions. It can be seen that the LDMD 
can be accurately determined without proper motions if the sam- 
pling region is the very narrow cone. On the other hand, in the 
case of the narrow cone, neglecting the proper motion measure- 
ments causes overestimation of the LDMD. These results do not 
change even if small distance errors with A = 0.1 are applied. 
Accordingly, it can be said that if tracer stars are observed in a 
sufficiently narrow conical region, proper motion measurements 
are not necessary. However, it may be difficult or impossible for 
such a narrow cone to contain a sufficient number of tracer stars; 
our result in Figure 2] suggests that the required sample size is 
6,000 stars. 

Data of a tracer population used in G12 w as taken from 
a K star catalog of Kuiike n & Gikno re fl989a) calibrated by 
[a/Fe] distributions from Hipparcos (Kotonevaet al. 2002) and 
Sloan Extension for Galactic Understanding and Exploration 
(SEGUE). Although their data do not seem to have proper mo- 
tion measurements, since their sampling region is as narrow a 
cone as our very narrow cone, the high LDMD determined by 
G12 cannot be attributed to ignoring proper motions. 

6.2.4. Line-of-sight velocity errors 

Since we cannot obtain LOSVs of stars from astrometry alone 
and need independent spectroscopic observations, observational 
errors on LOSVs have no relation to the distance and proper mo- 
tion errors. As well as the functional forms of the astrometric 
errors we introduced with Eq. (l2ll and (l22t . here we formalize 
the LOSV error as 

^''"""^(li^) (2^^ 

where C and /3 are parameters. C corresponds to the LOSV un- 
certainties at d — I kpc in the unit of km s^'. If quality of the 
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spectroscopic observations is good enough, the LOSV precisions 
would be Hmited to wavelength resolutions and/or calibration er- 
rors of the observations: precision floors. In this case, the LOSV 
precisions can be assumed to be independent from distances to 
stars, which correspond to j3 = 0. Otherwise, we assume that 
the LOSV errors depend on signal-to-noise ratios of stars. If 
the signal-to-noise ratios are dominated by photon statistics, the 
LOSV errors can be ideally assumed to vary with stellar flux as 
£iosv °^ 1 / v(/ °^ ^- This case corresponds to )3 = L However, 
the L OSV errors in ac tual obser vations are h ighly complicated 
(e.g.. lKatz etai]|2004l :[de Bmiin e et alj|2005h . We therefore ex- 
amine cases in which j3 = 0-3 is adopted; in actual LOS V mea - 
surements, j3 seems to be in the range of 2-3 (e.g., Prustij l2012l) . 
However, again it should be noted that this formalization is sim- 
plistic. Even if it is the case of j3 > 1, the uncertainties are limited 
to the precision floors at li ~ 0. Besides, the uncertainties depend 
on stellar populations. Because the LOSV measurements are in- 
dependent from astrometric observations, we assume that there 
are no astrometric errors here: A = Z? = 0. In this case, the tracer 
density is not affected by the observational errors. 

Figures[TO]and[TT]show our results of adopting the various C, 
j3 and sampling regions. In the cases of j3 = 0, the conical sam- 
pling regions overestimate the LDMD when C > 3.0 kms^', 
and the cylinders overestimate when C > 5.0 km s^^ in some 
runs. These values of C can be considered as the required pre- 
cision floors of the LOSV measurements. Where j3 = 1, all of 
the sampling regions overestimate the LDMD when the eiTors 
are > 7.5 km s^^ at c/ = 1 kpc in some runs. Where j3 = 2 and 
3 behavior of derived LDMDs is interesting; the LDMDs are 
underestimated. Where = 2, the LOSV precision required to 
prevent significant underestimation is approximately 20 km s^' 
aid = \ kpc for the conical sampling regions and 30 km s^^ for 
the cylinders. Where j3 = 3, the precision required is approxi- 
mately 10 km s^' at = 1 kpc for the conical sampling regions 
and 15 km s^' for the cylinders. 

Figure [H] shows changes of a°}^^^{z) by the LOSV errors. 
The behavior of the velocity dispersion change is complicated. 
When j3 = 0, the very narrow cone and the wide cylinder are the 
most and the least susceptible to the LOSV eiTors, respectively; 
however, when j3 > 2, the very narrow cone and the wide cylin- 
der are the least and the most susceptible, respectively. Figure 
[T3] shows tracer density fall-offs deduced from Eq. (|9]l under the 
LOSV errors. All parameters of p,(0), 0^j{Q) and Pdm are set to 
the true values. This figure explains why the LOSV eiTors un- 
derestimate or overestimate the LDMD depending on j3 . Where 
j3 = and 1, the MA method tends to predict thicker distribu- 
tions of the tracer than in the case of excluding the LOSV errors; 
on the other hand, in the cases of /3 = 2 and 3, the MA method 
predicts thinner distributions. Generally, higher LDMDs make 
the tracer distribution thinner, while lower LDMDs yield thicker 
tracer distribution. Therefore, in the cases of j3 = and 1, high 
LDMDs are required to match with the observed density fall- 
offs; in the cases of j3 = 2 and 3, low LDMDs are preferred. 

7. Discussion 

7. 1 . Required precisions and comparison witii tiie previous 
studies 

In Sect. 16.21 we generated mock data of tracer stars and applied 
the MA method to them. Sect. 16.2.11 is aimed at estimating the 
required sample size of the tracer In both cases of the conical 
and cylindrical sampling regions, we find that the sample size re- 
quired to determine the LDMD with accuracy is approximately 



6,000 stars. Gl 1 have used three tracer populations in their study, 
and their sample sizes are 139 K giants, 2026 A stars and 3080 F 
stars, respectively. Although the total number of their tracer stars 
is comparable to the required sample size we estimated, each of 
them is insufficient in number Moreover, although the shortest 
sample of the K stars extends from z = 0.2 kpc to 0.7 kpc, the 
other two samples only cover the range of z < 0.2 kpc. Therefore, 
we suspect that the samples of Gil may lack sufficient tracer 
stars in high z regions. On the other hand, G12 used a single 
tracer population which consists of 2016 K stars used for the 
density fall-off and 580 K stars used for the velocity disper- 
sion profile. The sample size of G12 is much shorter than the 
required sample size we estimated. Therefore, we also suspect 
that the LDMD of G 12 may be imprecis ely determined due to 
statistical uncertaintv. ' Zhang et al.l (1201 2h more recently derived 
the LDMD using their own me thod; their me asurement of the 
LDMD was lower than that of Garb ari et al.l by a factor of 3 
or 4. They used three tracer populations containing 3672, 1416 
and 20 01 K dwarfs. As well as Gl 1, each tracer of IZhang et al.l 
is insufficient in number although the total sample size 
is larger than the required sample size. However, the method of 
IZhang et"aD is different from the MA method; therefore, we can- 
not simply compare their sample size with our result. 

In Sect. 16.2.21 we apply the astrometric parallax measure- 
ment errors to the tracer stars. The distance errors deform the 
profiles of density fall-off and vertical velocity dispersion. Our 
result shows such faulty profiles can overestimate the LDMD. 
We estimate that the required precision of the parallax measure- 
ments is approximately ^0.1 mas at li = 1 kpc in almost all 
cases. This required precision has not been achieved in present 
astrometric observations yet: Standard eiTors of Hipparcos ob- 
servations are in the range of = 0.7-0.9 mas even for stars 
brighter than 9th magnitude (mag) in a catalog published in 1 997 
(Ferryman et al. 1997; Ferryman 2012)Q Since Gil have used 
the Hipparcos data for density profiles of their tracers, their 
tracer data are less precise than the required value. Although 
Gil did not use observed velocity data but introduced paramet- 
ric runs for 0?\^^^{z), as we showed in Sect. 16.2.21 the LDMD 
can be overestimated solely by deformation of the tracer den- 
sity profile caused by the distance errors. Therefore, we in- 
fer that the LDMD determined by Gil could be overestimated 
due to the imprecise astrometric measurements. Since G12 
and Zhang et al. (2012?) used photometric distances provided 
by Kuiiken & Gilmora (Il989ah and Sloan Digital Sky Survey 
(SDSSySEGUE, respectively, we cannot perform fair compari- 
son with them. 

It should be noted that our results may depend on details 
of our settings. For example, we arbitrarily set the scale height 
of our tracer to /inac = 200 pc. Besides, Zniin — 0.2 kpc and 
Zmax = 1.2 kpc are also arbitrary. In Gil, their samples cover 
up to the height of Zmax — 0.75 kpc. It is expected that such a 
short sample can avoid high z regions where the observational 
errors become large. Figure [T4]shows determined LDMDs in the 
cases of setting Zmax — 0.75 kpc. Here, the parallax measure- 
ment eiTors are taken into account. Although the overestimation 
of the LDMD is slightly mitigated in the case of the naiTow cone 
sampling region, significantly high LDMDs are predicted when 
A > 0.5. Therefore, we still expect that the use of the Hipparcos 
sample causes overestimation even if tracer stars in the high z 



Although Ivan Leeuwen & Fantind ( l2005h significantly improved 
the accuracy of the Hipparcos data, (7 ^ < 0.3 mas can only be ach ieved 
for stars brighter than 6 mag (see also Ivan Leeuwenll2005L 12007 al lbh. 
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regions are excluded. Our results conclude that the required par- 
allax precisions are an; = 0.1-0.3 mas. 

In Sect. 16.2.31 we estimate the required precisions of proper 
motion measurements, which seem larger than 3 masyr^' at 
d = 1 kpc. In the cases of the cylindrical sampling regions, 
however, the required precision becomes severe: approximately 
(7^ < 1 mas yr^ ' at = 1 kpc for the narrow cylinder and G/i < 
0.3 m as yr~^ for the wide cylinder, respectively. Zhang et"aD 
(1201 2h sampled their tracers in a wide cylinder the radius of 
which is 1 kpc. Their proper motion data taken from the SDSS 
catalog are as precise as ~ 3 mas yr~' for stars brighter than 
18 mag in r-band (e.g.. lBond et alJl2010h . Their proper motion 
precisions would be insufficient for the wide cylinder sampling 
region ;/ the MA method is used. However, since Zhang et al. 
(1201 2t) determined the LDMD with their own method, we can- 
not conclude that their result is overestimating. 

In Sect. 16.2.41 we estimate the required precisions of LOS V 
measurements and find that the LOSV eiTors lead to either 
overestimation or underestimation of the LDMD, which is con- 
tingent on distance-dependence of the errors. When we apply 
uniform errors: fi — Q, the required precision is found to be 
< 3.0 kms^^ for the cylindrical sampling regions and 5.0- 
7.5 kms^^ for the cylinder sampling regions. The required 
precision in this case coiTesponds to required precision floors 
in observations. Where the errors only weakly depend on dis- 
tances, the LOSV errors can cause overestimation of the LDMD; 
whereas, in the cases of strong distance-dependence of the er- 
rors, the LO SV errors can cause underestimation. G12 used 
the sample of iKuiiken & GilmorS ([T989a) . the LOSV errors of 
which are limited to < 10 km s^' at J = 1 kpc. Therefore, their 
LOSV data would not cause underestimation or overestimation 
of the LDMD as long as )3 > 2 and their precision floors are 
sufficiently low. 

To summarize our results, the parallax and proper mo- 
tion uncertainties we examined can cause overestimation of the 
LDMD in calculation using the MA method. Gil have used the 
Hipparcos data, the parallax measurements of which are not pre- 
cise enough to determine the LDMD accurately. On the other 
hand, G12 used the photometric distances. Therefore, we cannot 
fairly evaluate the reliability of their result. However, since G12 
also resulted in the high LDMD which is consistent with Gl 1, it 
is expected that G12 would also be overestimating. The LOSV 
data of Gl 1 and G12 seem to be sufficiently precise. It is now in- 
teresting to return to Figure[T] In the right panel, the PDF of Pdm- 
Es inclines from the upper left to the lower right; the same trend 
can be seen in Figure 8 of G 12. This means that pdm and Es have 
a coiTelation: When Pdm is overestimated, Ej would be underes- 
timated. Figure[T5]shows the PDFs in some cases where LDMDs 
are ill-determined. As we expected. Eg is underestimated when 
Pdm is overestimated and vice versa. In the results of G12, a bary- 
onic surface density they determined is Es ~ 45.5^5 ^ M© pc^^, 
and an observed surface density they applied in their calcula- 
tion is Ef' = 49.4 ±4.6 Mq pc^^ (Flynn et al. 2006). Although 
the median value of the calculation is in the range of observa- 
tional uncertainty, there could be small underestimation of Eg. 
This may be a hint concerning the overestimation of the LDMD 
in G12. However, we emphasise our result that the MA method 
itself does not indicate any systematic biases and is capable of 
determining the LDMD with accuracy if sample size and obser- 
vational precisions are sufficient (Sect. l6Tt . 

Our results, however, do not necessarily mean that the rel- 
atively low LDMDs suggested by the other previous studies 
are approved. Quality of the data used in Gil and G12 are 



not considerably inferior to the others. Although G i l and G12 
demonstrated that methods of Hol mberg & Flvniil (l2000l) and 
Kuijken & Gilmore ( 1989a b e) have systematic bia ses to under- 
estima te the LDMD, the other methods (e.g., the IZhang et"al] 
(12012!) method) and accuracies of their data should also be as- 
sessed carefully. F or example, photometric distances may have 
systematic errors; iBovv et aH (1201 2dh have shown that stellar 
distances of G dwarfs determined with an isochrone of A n et all 
(|2009|) are n early 10 % larger than those determined with an 
isochrone of l lvezic et al.l (l2008h . Systematically large distances 
of tracer stars result in large scale heights of the distributions. 
Such thick tracers can underestimate the LDMD. 



7.2. Validity of our galaxy modeling 

We have to mention, however, that our analytical model to gen- 
erate the mock data would be somewhat simplistic. For example, 
we assumed the galaxy model consisting of the 15 visible mat- 
ter components. Such mode ling may not n ecessarily well repre- 
sent the real Galaxy. Bow et al.l (l201 2b"c.d) have recently dis- 
cussed that the Galactic disk is composed of a lot of continuous 
"mono-abundance sub-populations" and that there are no dis- 
tin ct thin/thick disks. If this is the case, the stellar classification 
bv lFlvnn et al.l (l2006l Table[3]in this paper) may not be suitable. 
In addition, as an observational difficulty, even if data are accu- 
rate and precise enough, it may be laborious to pick out kinemat- 
ically homogeneous tracer stars and to verify the homogeneity. 
Contaminated tracer samples can lead to ill-determined LDMDs 
(e.g., Kuijken & Gilmore 1989a b). This also needs reliable stel- 
lar population synthesis theory and spectroscopic observations. 

The tracer in our study was modelled to have the uniform 
density and velocity distributions in a z-plane. Actually, stellar 
density and velocity dispersions in the Galaxy vary in the radial 
direct i on with scale radii of ^3.0 kpc (e.g., Lewis & Freemai^ 
119891; iBovv et all l2012dl) PI Therefore, if tracer-sampling re- 
gions are too wide, the kinematic homogeneity for the tracer can 
be broken in actual observations (see Gl 1). 

If our modeling can represent the real Galaxy well, our 
mock observations may be used to calibrate erroneously deter- 
mined LDMDs. An advantageous point of astrometric obser- 
vations is that observational uncertainties are already known. 
Therefore, even if sufficiently precise data are not available, our 
method can estimate magnitudes of under/overestimation of the 
ill-determined LDMDs by applying the known astrometric er- 
rors to the mock data. 



7.3. Toward the future astrometric observations 

Our results have shown that the MA method requires a high de- 
gree of astrometric precision to accurately determine the LDMD, 
which still cannot be achieved by present astrometric observa- 
tions. The near-future astrometric satellite Gaia, however, will 
significantly exceed Hipparcos. Gaia is expected to observe the 
complete sample of all stars brighter than 20 mag with end-of- 
mission parallax precisions of ~ 0.01 mas at V = 10 mag (V de- 
notes Johnson V magnitude), 0.01-0.03 mas at V = 15 mag and 
up to 0.1-0.35 mas at V = 20 mag (e.g..lde Bruiine et al.ll200l 



IBailer-Jonesl2009l; iJordi et al.l20Tot lPrustil20 1 2h ; proper motion 
precisions in the unit of mas yr^' are comparable to the parallax 
precisions. A frequent tracer population is K dwarfs, t he abso- 
lute magnitude of which are V ^ 7 mag (iBinnev & Merrifieldl 



iBovv et aLl ( l2012cl) have observationally shown that the scale radii 
of vertical velocity dispersions are as long as ~ 7.1 kpc. 
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Il998l) . This means that K stars at c/ = 1 kpc have the bright- 
ness of y ^ 17 mag. Therefore, these stars are expected to be 
observed with the precisions of < 0.1 mas. Thus Gaia will be 
able to achieve the required astrometric precisions we estimated. 
Furthermore, Gaia is designed to measure LOSVs of the stars 
simultaneously. The LOSV precisions are 10-1 7 km s~^ at V = 
17 m ag with precision flo ors of ^ 1 km s~' (Ide Bruiine et al.l 
[2005; Katzetal. 2004; Jo rdi et all 1201 Oh . distance-dependence 
of which seems to be j3 = 2-3 (from comparison with Figure 7 
of lPrustill2012h . From comparing with Figure [TT] Gaia seems 
to be able to determine the LDMD within the ranges of 90 % 
confidence levels although there may be slight underestimation 
if j3 = 3 and C = 15. Accordingly, Gaia is expected to enable 
determination of the LDMD using the MA method although it 
may be better to replace the LOSV data with more precise mea- 
surements. 

In addition to Gaia, Nano-JASMINE and JMAPS are also 
individually planned to produce all-sky astrometric catalogs. By 
combining their data with the Hipparcos and Tycho positional 
data from 1991, proper motions with 0.1 mas yr^' are expected 
to be achievable for bright stars (Michaliket al. 2012). Such a 
high precision astrometry would improve not only the LDMD 
determination but also our understanding of various aspects of 
Galactic dynamics. 

8. Conclusions 

We scrutinized the MA method devised by iGarbari et al.l and 
carefully assessed their LDMD determinations. We created 
mock observational data and applied the MA method. As re- 
sults, we found that the MA method is capable of determining 
the LDMD with accuracy if the sample size of a tracer and ob- 
servational precisions are sufficient. We found that the sample 
size must be larger than approximately 6,000 stars. Astrometric 
errors, however, can cause overestimation on the LDMD deter- 
mination, and we estimated that the required precisions are ap- 
proximately 0.1-0.3 mas for parallax measurements. Proper mo- 
tion precisions do not seem to be as important as the parallax 
precisions. In addition, LOSV errors can cause either overesti- 
mation or underestimation of the LDMD determination: If the 
errors weakly or strongly depend on distance, the LDMD can be 
overestimated or underestimated, respectively. Although the re- 
quired precisions of LOSV measurements are also contingent on 
the distance-dependence of the errors (see Sect. I6.2.4l and Figure 
[To] andOfor details). 

From our results, we infer that Gil who used Hipparcos 
data would be overestimating the LDMD due to the insuffi- 
cient astrometric precisions. Moreover, although the LDMD de- 
termined by G12 who used photometric distances is consistent 
with the result of Gl 1, the sample size of G12 would not be large 
enough to precisely determine the LDMD. However, our results 
do not necessarily mean that other previous studies displaying 
a preference for relative ly low LDMDs (e.g., Zhang et al. 2012; 
iBovv & Tremainel2012l) are satisfactory. Careful scrutinizations 
for the other methods to determine the LDMD are also needed. 

The near-future astrometric satellite Gaia can be expected to 
measure parallaxes, proper motions and LOSVs with sufficient 
precisions. Our results indicate that Gaia will provide data that 
enables us to determine the LDMD using the MA method. 
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overestimate in z > 0.2 kpc. However, magnitudes of the over- 
estimation at Zmax are only ^ 0.5 % and ^ 1.5 % in the cases of 
the isothermal and the quadratic dispersions, respectively. From 
these results, we can see that the isothermality for the deduced 
model hardly affects the determination of the LDMD. In addi- 
tion, the runs of velocity dispersions of the model components 
do not seem to be important. Gil have also come to the same 
conclusion and they mentioned as follows: Although the veloc- 
ity dispersion profile of the tracer directly affects the tracer den- 
sity fall-off and has significant impact on the result, uncertainties 
in dispersion profiles of the model components are, by contrast, 
marginalized out when we calculate p^m and Ps(0), since they 
appear only in Eq. (O through Eq. (|3]l. From the above results, 
it is verified that assuming isothermality for the model compo- 
nents does not lead to erroneous determinations of the LDMD. 
Although the quadratic function used in Gil is not necessarily 
better than the isothermal approximation, the impact of the dif- 
ference between them seems trivial. 



Appendix A: Validity of the isotliermality for the 
deduced model 

Deduced galaxy models of Gil and G12 are slightly different 
from each other Gil has used parametrized runs of velocity dis- 
pe rsions assurning a behavior similar to the observational fitting 
bv lBondet al.1 ([2010). G12, on the other hand, assumed that all 
components are isothermal. Here, we discuss the impact of these 
velocity dispersion profiles on the calculation of the MA method. 
In addition, we verify the appropriateness of the isothermal as- 
sumption used in G12 and this study. 

The dispersion runs introduced in Gil is described by a 
quadratic form: 

(T?,(z) = (t2,(0)(1+cz2), (A.l) 

where the constant of c is chosen so that thi s function satisfies 
the observational fitting o f iBond et al.l(l2010h : 

CT,,KZmax) = (Tz,,(0)+4f^j . (A.2) 

Figure lA.ll shows comparison between the quadratic functions 
and the analytical solutions of the Jeans equation for some model 
components. Although some components seem to be fitted well 
by the quadratic functions, some others are fitted poorly. This 
would be due to the fact that the observations of iBond et alJ 
(l20Toh are for blue disk stars and/or our assumed model is ar- 
tificial. 

We set all parameters of Pj (0), o;,, (0) and p^m to the true val- 
ues and compute surface density profiles of the deduced models, 
(47rG)^'d<I>/dz, following the procedure of Sect. 12. II with the 
quadratic functions and the isothermality. Figure |A!2] shows the 
results. In the top panel, we compare the results with the as- 
sumed density profile of Eq. (fTTl i. Since all of them are almost 
consistent, it can be said that both the quadratic functions and 
the isothermality can reproduce well the assumed galactic po- 
tential. The bottom panel shows fractional differences between 
the results and the assumed model. In both cases, the deduced 
models underestimate the surface density in z < 0.2 kpc and 
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Fig. 9. LDMDs determined in the cases of ignoring the proper 
motion measurement. Here we set ^l0=O for all tracer stars. 
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Fig. 8. Influence of distance and proper motion errors on vertical 
velocities of the tracer stars in three example cases. In each three 
panels, the top left panel compares true and those affected 
by the distance errors alone. The top right panel compares the 
affected by the distance errors and those affected by both the dis- 
tance and the proper motion errors. The bottom left panel shows 
the true and those affected by both errors. 
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Fig, 12. Changes of vertical velocity dispersion profiles of the 
tracer by the LOSV errors. Here we used mock data of one mil- 
lion stars. 
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Fig. 13. Tracer density fall-offs deduced by the MA method. 

Here we set all parameters to the true values and apply the LOSV 
errors only. Because distance errors are not taken into account, 
the deduced profile in the case of no errors corresponds to the 
observed density fall-off. The mock data of one milhon stars are 
used. 
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Fig. 15. PDFs of LDMD, local and surface densities of the vis- 
ible matters. The red and green lines are the same as in Figure 
[T] The top panel shows the result in the case of the narrow cone 
sampling region with A = 0.3, B = C = Q, a = Q and the sample 
ID = 1. This case overestimates the LDMD. The middle panel 
shows the result in the case of the wide cylinder sampling region 
with A = 0.1, B = 2.0 C = 0, a = and the sample ID = 1. This 
case overestimates the LDMD. The bottom panel from the top 
shows the result in the case of the narrow cylinder sampling re- 
gion with A=B = Q,C = 30.0, j3 = 3 and the sample ID = 1. 
This case underestimates the LDMD. 
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Fig. A.l. Velocity dispersion profiles of some model components 
in Table [T] The solid and dashed lines indicate analytical solu- 
tions computed by Eq. (fTSl i and fittings by the quadratic function, 
respectively. 
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Fig. 5. Influence of the astrometric distance errors on determining the LDMD. The left and right panels indicate the results of 
adopting a = and a — 1, respectively. In each case, we conduct the same computation with different tracer samples generated in 
tile same conditions (ID = 1-6). The horizontal dotted line indicates the true LDMD. In the case of the wide cylinder with a = 1 
(tiie right bottom panel), we additionally show the result of adopting A = 0.05. 
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Fig. 7. Influence of the proper motion errors on determining the LDMD. The left and right panels indicate the results of adopting 
a = and a = 1 , respectively. In all cases, we set the magnitude of distance error to be A = 0. 1 . In each case, we conduct the same 
computations with different tracer samples generated from the same conditions (ID =1-6). The horizontal dotted line indicates 
the true LDMD. In the case of the wide cylinder with a = 1 (the right bottom panel), we omit the result of adopting B = 3.0 since 
derived LDMDs become extremely high. 
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Fig. 10. Influence of the LOS V errors on determining LDMD. The left and right panels indicate the results of adopting j3 = and 
/3 = 1, respectively. In all cases, we assume that there are no distance and proper motion errors: A = B = 0. In each case, we conduct 
the same computations with different tracer samples generated from the same conditions (ID =1-6). The horizontal dotted Une 
indicates the true LDMD. 
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Fig. 11. Same as Figure[TO]but adopting j3 = 2 (left panels) and 15—3 (right panels). 
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